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We will present the numerical evaluation of the hopping and mag- 
netic exchange integrals for a nearest-neighbor t — J model of the 
quarter-filled a'NaV^Os compound. The effective integrals are ob- 
tained from valence-spectroscopy ab tnitio calculations of embedded 
crystal fragments (two VO5 pyramids in the different geometries cor- 
responding to the desired parameters). We are using a large config- 
urations interaction (CI) method, where the CI space is specifically 
optimized to obtain accurate energy differences. We show that the 
a'NaV205 system can be seen as a two-dimensional asymmetric trian- 
gular Heisenberg lattice where the effective sites represent delocalized 
V — O — V rung entities supporting the magnetic electrons. 



I. INTRODUCTION 

Spin-Peierls (SP) transitions have been long known in organic com- 
pound Jil but the observation, in 1993, of a SP transition in the inorganic 
CuGeOsP {Tsp = 14:K) renewed the interest of the physical commu- 
nity for this phenomenum. In 1996 the a'NaV20M inorganic system 
was found to present a second order SP transitiona at Tsp = SAK, the 
highest Tsp temperature so far known. Indeed a rapid drop of the mag- 
netic susceptibility, the opening of a spin gap! (A = 9.8meV) as well as 
the magneto-distortion of the lattice was observed in this compound. 
This compound has attracted special attention since, in addition of 
its very high Tsp, the value of the j^^^^ = Q-4:4: ratio does not agree 
with the BCS theory predicted value of 3.53, the Tsp dependence with 
the magnetic field is much weakeiil than the theoretical predictions! 
and the temperature dependence of the thermal conductivity is much 
larger than in CuGeO^p. 

The crystal is formed by layers of VO5 square-pyramids stacked along 
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the c axis. The oxygen atoms of the pyramid basis form a quasi-square 
planar lattice along the a and h directions. The pyramids are alterna- 
tively pointing on top or below these planes. Periodical vacancies of 
the VO top of the pyramids are replaced by Na^ ions forming chains 
parallel to the b axis (see fig. |T]). 

For a long time, the high temperature crystal structure was assumed 
to be non-centrosymmetric (in P2imn group )§, and consequently the 
electronic structure was assumed to be an alternation of half filled V'^'^ 
magnetic chains and of non-magnetic V^~^ chains along the b axis. 
The system was then supposed to be constituted of isolated Heisen- 
berg chains and the SP transition was naturally understood in this 
description by the one- dimensional dimerization of the V'^'^ magnetic 
chains. 

However the existence of a pseudo-inversion center in this compound 
led a couple of groups0 to re-investigate the crystallographical data us- 
ing more recent experimental methods. They found that a'NaV20:i was 
centrosymmetric (in the Pmmn space group) destroying the grounding 
for the assumed differentiation between V'^'^ and V^'^ chains. a'NaV205 
should then be seen as a quarter-filled system. 

It is widely accepted, after GalyS, that the magnetic interactions take 
only place within the layers, and thus the pertinent model for a'NaV205 
is a quarter-filled t — J model Hamiltonian in the (a, b) plane, based on 
the V magnetic 3d orbitals, where the interactions are limited to the 
nearest neighbors (NN) magnetic sites. 

<i,j> o" <i,j> ^ / 



— * 

where the sum over < i,j > runs over NN magnetic sites. Si is the 
local spin operator on site i, (resp. Cj^) are the usual creation (resp. 
annihilation) operators of an electron of spin a on site i, rij is the 
number operator on site z, Jij is the magnetic exchange integral and Uj 
the hopping integral of a magnetic electron between sites i and j. 

In the present system one has three different types of NN interactions 
(see fig. 0) : the interactions along the ladders (in the h direction) 
denoted as t\\ and Jy, the interactions along ladders rungs denoted as 
t_L and J_L (in the a direction) and the interactions between the ladders 
denoted as t' and J' . According to the order of magnitude of these 
parameters, the system electronic structure can be very different, going 
from a one dimensional behavior with essentially non interacting two- 
legs ladders (|t'|, | J'| <^ 1), or zig-zag frustrated ladders (|t±|, | J±| <^ 
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1) to a two-dimensional system. It is therefore of crucial importance 
to have a reliable evaluation of the relative amplitude of the different 
parameters. 

The aim of this work is to evaluate, using ab initio quantum chemi- 
cal methods, the effective hopping and exchange integrals between NN 
vanadium atoms. 



II. METHOD 



The exchange and hopping effective integrals are essentially local pa- 
rameters, therefore they can be accurately evaluated using fragment 
spectroscopy, provided that (i) the fragment includes all the crystal 
short range effects, that is the local environment of the magnetic atoms, 
and (ii) the crystal long range effects are treated through an appropriate 
bath. 

We performed excitation energies calculations on bi-pyramidal frag- 
ments (including the two NN magnetic atoms as well as all their sur- 
rounding oxygens) in a bath reproducing the main effects of the rest 
of the crystal (Madelung potential and exclusion effects). The frag- 
ments used for the t\\ (resp. J\\) and t± (resp. J±) calculations are built 
from two VO5 pyramids sharing a corner (see fig. ^ a and b) while the 
fragment used for the t' (resp. J') calculations is constituted of two 
VO5 pyramids sharing an edge (see fig. ^ c). Two sets of calculations 
are performed for each fragment. The first one involving two unpaired 
electrons in the two VO5 pyramids allows us to obtain the effective ex- 
change parameters (including all direct and super-exchange processes) 
from the singlet-triplet first excitation energy. 

J = Es-Et (2) 



The second set of calculations, involving only one unpaired magnetic 
electron in the two VO^ pyramids, allows us to obtain the effective 
hopping parameters from the first doublet-doublet excitation energy. 

t = E£±^^ (3) 



where the (resp. D_) doublet is of the same symmetry than the 
3di + 3d2 (resp. 3di — 3^2) delocalized orbital on the two magnetic 
centers. One notes that according to the relative energies of the 
and states, the sign of the hopping integral can change. 
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A. Modeling the rest of the crystal 



As state above a correct fragment calculation should take into account 
the main effects of the rest of the crystal. In our ionic, strongly corre- 
lated system these effects are limited to the short range Pauli exclusion 
effects and the long range Madelung potential. 

First, the Madelung potential is reproduced by a set of positive and 
negative charges corresponding to the cations and anions. We used a 
cutoff threshold associated with a ISA radius sphere, centered on the 
bi-pyramid fragment. The border of the sphere is designed so that 
to preserve VO5 pyramids, the chemical meaningful entities. Border 
charges are adapted by an Evjen procedure. 

The a'NaV205 system presents 4 different types of atoms : the fully 
ionized sodium cations, the mixed valence vanadium atoms and two 
types of oxygen atoms, the in-plane oxygens and the apical oxygens. 
These two types of oxygen play a very different role in the chemistry 
of the compound. Indeed, ab-initio calculations on the bi-pyramidal 
systems, both at the Hartree-Fock (HF) and correlated level, as well 
HF calculations on the whole crystal give the same result, the apical 
oxygen is bonded to the vanadium atom by a slightly polarized double 
bond. The Miilliken population analysis are similar in all calculations 
with a charge close to — 2e for the in-plane oxygens while the apical ones 
present only a small charge of about — 0.5e. In fact the VO^ pyramids 
should not be seen as such, but rather as {V = 0)^+°''^+ molecules on 
an oxygen double-anions square lattice. This picture is confirmed by 
Raman experimentsEl. Indeed, two distinct bands are observed, one for 
the stretch of the vanadium/apical-oxygen bond located at 970 cm~^, 
the other for the group of vibrations of the vanadium against the oxygen 
plane which is located in the 400-480 cm~^ range. 

The embedding charges are then taken as follow : 

• the sodium atoms are represented by +1.0e charges, 

• the oxygen atoms of the quasi-square lattice forming the basis of 
the pyramids are represented by — 2.0e charges, 

• the apical oxygens are represented by — 0.5e charges 

• and the vanadium atoms are represented by their average valence 
of +3.0e. 

Second, the exclusion effects are taken into account through total 
charges pseudo-potentials. The embedded fragments are highly nega- 



4 



tively charged (— 8e or — 9e) and their electrons would like to expand 
out of the fragments volume, thus the contention effect of the rest of the 
crystal needs to be modeled by an exclusion potential. The negative 
charge of the oxygen anions, even modeled as point charges, is sufficient 
to insure the exclusion effect on the fragment electrons ; however the 
positive charge of the cations attracts the fragment electrons and it is 



necessary to explicitly treat the exclusion, 
the cations with total-ions pseudo-potential^'^^ 



lis is done by modeling 
using the Durand and 



Barthelat large core effective potentialsB, instead of only point charges 

The good behavior of the above embedding has been checked against 
HF crystal calculations. Figure § shows the Projected Density of States 
(PDOSS) on the magnetic orbitals of the vanadium atoms (1) and on 
the p orbital of the bridging oxygen of the rungs (2). The PDOSS 
is compared with its equivalent in the cluster calculations, the square 
of the cluster orbitals projection onto the atomic orbitals, plotted as 
a function of the cluster orbital energies. On can see that while the 
above embedding yields a favorable comparison of the cluster orbitals 
with the crystal PDOSS, the embedding defined with only one kind of 
oxygen atoms and the formal charges (0~^, V^"*"^'^, Na~^) do not properly 
place the vanadium d orbitals, in particular the gap at the Fermi level 
is overestimated. In order to emphasize the importance of a proper 
embedding, the isolated cluster orbitals are also reported. The result 
goes without comment. 



B. Computing accurate excitation energies 



Ab initio quantum chemical methods are powerful tools to obtain reli- 
able excitation energies as well as good representation of the associated 
ground or excited states wave functions. As will be seen later a precise 
analysis of the wave function informations allows a deep understand- 
ing of the different mechanisms involved in the effective exchange and 
hopping processes. 

Let us first consider the physics of the excitations we are trying to 
compute. One sees easily that the fragment orbitals can be divided into 
three subsets according to their role in the many-body wave functions 
of the two states we are looking for. 

The occupied orbitals are the orbitals that always remain doubly- 
occupied in the many-body states. 

The active orbitals are the ones which occupation number can 
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change in the different Slater determinants participating to the 
many-body states. In our case these orbitals are not only the two 
3d magnetic orbitals (from now on we will refer to them as di and 
d2 according to the V atom to which they belong) located on the 
two vanadium atoms but also the 2p orbitals (we will refer to them 
as p) of the bridging oxygen atom(s) that can participate to the 
effective exchange through super-exchange mechanism and to the 
hopping process by a through bridge delocalization of the hole. 
Symmetry considerations reduce the number of these bridging or- 
bitals to two 2p orbitals (one per irrep. involved). One should 
note that these active orbitals correspond to the valence orbitals 
that would support a three-bands Hubbard representation of the 
system — where both the V and the O atoms of the pyramids 
bases are considered. 

The virtual orbitals are the orbitals that always remain empty in 
the many-body states. 

A good "zeroth-order" description of the seeked states will then be 
provided by the eigensolutions of the exact Hamiltonian in the Com- 
plete Active Space (CAS) ; the CAS is defined as all possible con- 
figurations built from the occupation rules given above. In fact we 
simultaneously optimize the orbitals and the wave function coefficients 
in a Complete Active Space Self Consistent Field (CASSCF) proce- 
durelll. At this level of calculation, the polarization and correlation 
effects within the active space (static polarization and correlation) are 
variationally taken into account, while the other electrons are described 
within a mean field approximation. In the J calculations, six electrons 
— the two unpaired electrons plus the four electrons originating from 
the bridging p orbitals of the central oxygen(s) — are distributed in all 
possible manners into the four active orbitals. In t calculations, only 
five electrons are distributed into the four active orbitals. 

"Dynamical" polarization and correlation effects coming from exci- 
tations out of the active shell (originating from the occupied orbitals 
or ending in the virtual ones) are however crucial in order to obtain 
reliable results for excitation energies. The "dynamical" polarization 
effects come from the single-excitations on the CASSCF wave function, 
while the dominant contributions to "dynamical" correlation effects 
come from the double-excitations on the CASSCF wave function. We 
will take these effects into account in a two steps procedure, 
a) A self-consistent procedure is performed. It involves the diagonaliza- 
tion of the CAS + all single excitations space as well as a new orbitals 
optimization at this level. This is the so-called Iterative Difference Ded- 
icated Configuration Interaction level 1 method (IDDCI1)0. 
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b) In a second step the double excitations on the CAS wave-function 
that participate to the energy difference (at the 2"'^ order of pertur- 
bation) are added to the CI space, which is diagonalized. This is the 
Difference Dedicated Configuration Interaction level 2 (DDCI2)tZl. 

The above method is powerful and reliable. It has been used with 
great success in a slightly different version (DDCI3) for the determina- 
tion of the t — J model parameters, for instance, on copper-oxideJll or 
ladder compounds such as SrCu20'^. The main differences between 
the procedure used in the present work and the works cited above are 
i) the optimization of the orbitals at the CAS + single excitations level 
(IDDCIl) instead of only at the CASSCF level ii) the inclusion of the 
dynamical polarization and correlation of the bridging oxygen p orbitals 
is done in a slightly different way. We choose to include the active or- 
bital(s) of the bridge in the CAS so that the dynamical polarization 
effects are taken into account in the DDCIl space and the dynami- 
cal correlation effects are taken into account in the DDCI2 space. De 
Graaf et a£l as well as Calzado et allll made the choice of a smaller 
CAS including only the two magnetic d orbitals and to extend the di- 
agonalization to the 2-holes 1-particle and 2-particles 1-hole CI space 
(DDCI3) in order to describe the same physics. The main advantages 
of our choice are the following 

• The optimization of the orbitals at a level including the dynamical 
polarization effects. 

• An equal treatment of the magnetic d orbitals and the active p 
orbitals of the bridge. This point can be important in the cases 
where the through bridge processes are large. 

• A much smaller size of the largest CI space to diagonalize. For 
instance in the calculation of the J±_ effective integral our choice 
yields a CI space of 592, 000 determinants while the DCCI3 space 
with the small CAS would give a CI of 2, 783, 000 determinants. 

In order to analyze the relative contributions of dynamical polariza- 
tion and correlation we performed an additional Configuration Interac- 
tion, namely the diagonalization of the CAS but defined on the set of 
orbitals optimized within the IDDCIl procedure (CASCI). 



C. Computational details 

All cluster calculations are performed using the Barandiaran basis 
setS with the recommended contraction: ?>s?>pAd for V atoms and 
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2s4:pld for O atoms. The effects of the inner electrons are modeled 
by core pseudo-potentialsEj, [ls^2s^2p^3s^] for V atoms and [Is^] for O 
atoms. 

The crystal HF calculations were performed using the CRYSTAL98 
packageBll. In order to be able to perform the infinite crystal calculation 
it has been necessary to reduce the basis set compared to the cluster 
calculations. We therefore used a basis setc3 of quality single-zeta for 
the core electrons, while the valence electrons were described using 
a double-zeta, plus one polarization function for the oxygen and two 
polarization functions for the vanadium. 



III. RESULTS AND DISCUSSION 

The computed hopping and exchange integrals are summarized in 
table |. 

The differences observed between the CASCI and the IDDCIl results 
point out the crucial importance of the inclusion of dynamical polariza- 
tion effects in order to obtain reliable results for effective exchange and 
hopping processes. In comparison the dynamical correlation effects are 
of much lesser importance since they modify the integrals amplitude by 
a factor of only 10% to 20%. 

The computed wave functions are reported in the Appendix. Their 
analysis shows the strong contributions of 0~ configurations on the 
bridging oxygen in the _L geometries and moderate ones in the || ge- 
ometries. These configurations are negligible in the prime geometry, 
where their weight is smaller than 1%. This is easily understood by the 
quasi-orthogonal character of the angles between the Vanadium atoms 
and the two bridging Oxygens, this quasi-orthogonality strongly hin- 
ders the overlap between the V magnetic orbitals and the p orbitals of 
the bridging oxygens. Finally it results in the ferromagnetic character 
of J' and the weakness of the t' hopping integral. 

In the rung or orthogonal geometry, the contribution of the 0~ config- 
urations is particularly large with weights of 15% in the triplet state, 
25% in the singlet state and weights as large as 49% {D_) and 53% 
in the two doublet states. This tremendous delocalization of the 
magnetic electron on the p orbital of the bridging oxygen is responsible 
for the accordingly large hopping integral t±, indeed the direct (by op- 
position to the through bridge) hopping mechanism contributes only by 
a few meV to the total amplitude. Similarly the strong 0~ contribu- 
tions in the singlet and the triplet states induce a large super-exchange 
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mechanism responsible for the strong anti-ferromagnetic character of 
J^. Different reasons can be pleaded to explain these very large delo- 
calization of the magnetic electrons on the bridging oxygen. The most 
important is the relative positions of the pyramids with respect to the 
bridging O atoms. The vanadium magnetic 3d orbitals are orthogonal 
to the V = O bonds and the angle of the V = O bond to the oxy- 
gen plane is closed, the corresponding tilting of the vanadium magnetic 
orbitals toward the oxygen plane enhances considerably the overlap be- 
tween these orbitals and the p orbitals of the bridging oxygen. Another 
important factor is the short distance (1.83A) between the vanadium 
atoms and the bridging oxygen in this geometry. 

The II geometry is slightly different since two p orbitals (of different 
symmetries) of the bridging oxygen are involved in the through bridge 
processes. Their consolidate weights in the doublet states are 20% in 
D_ and 26% in and they result in a moderate enhancement of 
t\\. In the singlet and triplet states however, the 0~ configurations 
do not contribute in a significant way to the wave functions since their 
consolidated weight is smaller than 1%. The result is a very weak super- 
exchange mechanism resulting in a slightly anti-ferromagnetic effective 
exchange integral (Jy ~ J±/60). 

The dominant interactions are by far those taking place on the lad- 
der rungs (t± = —538.2meV and J± = — 293.5meV), thus they should 
determinate the main representation of the electronic structure. Two 
different descriptions can be proposed depending which of the exchange 
or the delocalization processes shall dominate. Either two electrons 
are paired in a singlet, one rung out of two (favored by the large 
value of J±), or an unpaired electron is delocalized on each rung (fa- 
vored by the large hopping integral tj_). A simple energetic calcu- 
lation shows however that the second solution is much more favor- 
able — with E2/Nrungs = t± = —538.2meV — than the first one — 
Ei/Nrungs = 1/2J± = —lATmeV. On the basis of such a representa- 
tion, the a'NaV20^ system can be seen as a two dimensional triangular 
Heisenberg system (see fig. |]) where the effective magnetic sites are de- 
localized on the V — O — V runes, in agreement with the representation 
suggested by Horsch and MackEl. A wave-functions analysis of the two 
many-body states involved in the t±_ calculations confirm this result 
since the determinants involving an unpaired electron on the bridging 
oxygen are as probable as the determinants where the magnetic elec- 
tron is located on the vanadiums. Indeed their relative weights are 1.00 
in the Z)_ and 1.21 in the D+. The exchange integrals between these ef- 
fective sites are then anti- ferromagnetic in the h direction, J'^^^ = J\\/2, 
and ferromagnetic between the effective chains (representing the lad- 
ders), J' = J' /A. 
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IV. THE EXTENDED HUBBARD MODEL 



A. Coherence of the t — J parameters 



One of the first reflexes in front of quantitative evaluations of effective 
integrals in a t — J model is to check the following relation 



A quick calculation, using table | results, shows that the above relation 
is far from being verified in our case. Why? 

Equation ^ comes from a perturbative evaluation, J = —At'^/Ud, of 
the exchange integral from an underlying Hubbard or extended Hub- 
bard model, supposed to be the exact representation. We saw in the 
previous section that the contributions of the 0~ configurations are 
tremendously large in the rung geometry, a Hubbard model based only 
on the magnetic d orbitals of the vanadium atoms is therefore not realis- 
tic and any representation pretending to be the reference should include 
the p orbital of the bridging oxygen (at least in the rung geometry). 

We will use the following notations: subscripts p and d respectively 
refer to the p orbital of the bridging oxygen and the magnetic d orbital 
of the vanadium. 5 = — ep is the orbital energy difference. Up 
and Ufi are the on-site repulsion energies. Vpd is the vanadium-oxygen 
nearest neighbor coulombic repulsion. The hopping integrals require a 
little more care since they can be encountered in different situations 
according to the number and nature of the surrounding electrons: 



where h refers to the mono-electronic part (kinetic and nuclear attrac- 
tion) of the Hamiltonian and V = l/ri2 is the bi-electronic repulsion. 
The p and d orbitals being localized on nearest neighbor atoms, the 
{pp\V\pd) and {dp\V\dd) bi-electronic integrals are a priori of the same 
order of magnitude as the mono-electronic part of the hopping {p\h\d) ; 
thus one should consider 4 different types of hopping integrals accord- 
ing to the number and the nature of the spectator electrons on the bond 
where the hopping takes place. Let us define 




(4) 



{p\H\d) 
{dp\H\dd) 
{pp\H\pd) 
{dpp\H\dpd) 



{p\h\d) 

{p\h\d) + {dp\V\dd) 
{p\h\d) + {pp\V\pd) 
{p\h\d) + {pp\V\pd) + {dp\V\dd) 
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t = {p\H\d) 
= {dp\H\dd) 
= {pp\H\pd) 
tP'^ = {dpp\H\dpd) 



Within such a formahsm the leading perturbative term of the effective 
exchange integral between the i and j vanadium sites comes at the 
fourth order of perturbation 



J\ — — 4 r5 8 T — (5) 

= 4Ji,^ + 8 J2,x (6) 



where = 6 - Up + Ud- Vpd and A2 = 26 - Up + 2Ud - Wpd- 

The first term in the expression of J± comes from the configuration 
(^\didipp) + Ippdjdj)^ /y/2 and while second comes from \dididjdj). The 
second term is usually negligible in front of the first one, while that 
one is equal to —At\/Ud- However, due to the strong contributions 
of the magnetic electron dclocalization on the bridging oxygen, this is 
not presently true. The ratio between Ji _l and J2^^ can be evaluated 
from the singlet wave functions coefficients. Indeed the second order 
perturbative expression of the singlet wave function is given by 

\dippdj) - \dippdj) 
1^5.)= ^ + 

\didipdj) — \didipdj) — \dipdjpdj) + \dipdjpdj) 

: — 2 + 

^ \did,pp) +pp\d^dj) r^'uJ^j\ 



V2 



where 



V2tP' 



d 



, _ 2V2tHP'^ 

Thus (4Ji,^)/(8J2,±) = V2C2/2C^ ~ 1.11 (see Appendix). It comes 
4J| ~ -lU.SmeV. 
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In the ladder direction the 0~ configurations have neghgible weights 
in the singlet and triplet states, thus the effective exchange integral 
can be expressed as usual : J\\ = 2K\\ — AtyUd = ^K^ + 4Jij|, where 
K\\ is the direct exchange integral. The latter is usually considered as 
negligible, however the extremely small value of Jy imposes to explicitly 
take these effects into account in the present case. A direct evaluation 
yields K\\ = 2.0meV and the super-exchange contribution comes to be 
4Ji II = d.lmeV. One can now properly verify the coherence of our 
model hopping and exchange parameters : 

f^' = 18.7 ^ ^ = 17.0 

^From the above analysis it is clear that any method that would 
not explicitly take into account the delocalization on the bridging oxy- 
gen would obtain doubtful evaluations of the effective integrals in the 
present case. Indeed Smolinski et found hopping integrals in agree- 
ment with our results using the LDA-I-U approach, however their ex- 
change integrals are very different to ours. This discrepancy can be 
attributed to the fact that they do not take into account the delocal- 
ization of the magnetic electron on the bridging oxygen that accounts 
for half of J± value. 



B. The extended Hubbard Hamiltonian 

In addition to the t — J model, it is possible to obtain from our energy 
and wave function calculations, the parameters for the extended Hub- 
bard model described in the previous section ; that is a model based 
on the d magnetic orbitals of the vanadium atoms and the p orbitals of 
the bridging oxygen in the rung geometry. The parameter extraction 
will be done using a least square fit method where both the energy 
differences and the wave functions coefficients of the extended Hub- 
bard model are fitted on the computed ab initio wave functions and 
energies, according to the intermediate Hamiltonian theorycl. The dif- 
ferent parameters of the Hubbard Hamiltonian are therefore optimized 
so that the Hubbard secular equations of the singlet, triplet and the 
two doublets are verified with the computed wave functions (taken as 
the normalized projection onto the CAS of the DDCI2 wave functions) 
and the computed energy differences. The resulting parameters in the 
rung geometry are summarized in table 0. 

These values are in reasonable agreement with the vanadium on site 
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repulsion obtained from the ladder geometry 

C/rf = 4(t||)V(2i^||-J||)=6.8el^ 



V. CONCLUSION 

We computed the effective hopping and exchange integrals of a i — J 
Hamiltonian for the high temperature phase of a'NaV20^. The param- 
eters evaluations were performed using an embedded cluster approach 
devised to properly take into account the electrostatic and exclusion 
effects of the rest of the crystal as well as dynamical polarization and 
correlation effects within the cluster. The results yield dominant inter- 
actions to be the delocalization of the magnetic electrons on the ladder 
rungs. This incredibly large delocalization is mediated by the appro- 
priate p orbital of the bridging oxygen. In fact the magnetic electrons 
should not be seen as supported by the d orbitals of the vanadium atoms 
in a quarter-filled system, but by orbitals delocalized on the 3 atoms of 
the rungs, that is on the d orbitals of the two vanadium atoms and the 
p orbital of the bridging oxygen. In such a representation the system 
is no longer quarter-filled but half-filled. The remaining interactions 
between these delocalized magnetic electrons devise a two-dimensional 
triangular Hciscnbcrg system where the effective exchange in the lad- 
der direction is anti-ferromagnetic while it is ferromagnetic in the two 
other directions. 

In addition to the t — J parameters evaluation, a thorough analysis 
of the ab-initio variational wave functions allowed us to evaluate the 
relative amplitudes of the underlying three-bands extended Hubbard 
Hamiltonian. The negligible difference obtained between the oxygen 
on-site repulsion and the p^d orbital energy difference. Up— 5 < 0.03eV, 
yields a quasi-degeneracy between the configurations where the mag- 
netic electron is located on one of the d orbitals of the vanadium atoms 
and the p orbital of the bridging oxygen. It results in the delocalization 
of the magnetic electrons on the rungs. 
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APPENDIX 



The wave functions obtained at the DDCI2 level are developed onto 
wave functions including up to 592000 slater determinants (according 
to the state and the geometry), nevertheless their major contributions 
are within the CAS space. The CAS configurations of the computed 
wave functions are reported below, after relocalization of the optimized 
CAS orbitals. This was done by the mean of a procedure based on 
Boys' methodH. 

In the rung geometry 

\didipdj) - \didipdj) - \dipdjdj) + \dipdjdj) 
2 

0.u \dApP) + \pPd,d,) ^ ^ 
v2 

+small terms . . . 



0.38 



\didipdj) + \didipdj) — \dipdjdj) — \dipdjdj) 



+small terms . . . 



^^,) = 0.62'^-^^'^ + '^^'^^-^ + 
v2 

-2\dipdj) + \dipdj) + \dipdj) 
u.oy ;= h 

Q ^ Jdjdjdj) + \didjdj) 
V2 

Q ^ Jdjdip) - \pdjdj) 
+small terms . . . 
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v2 
V2 

QQQ ldjdjdj) - \didjdj) 
Q ^^ Idjdip) + \pdjdj) 
+small terms . . . 



In the ladder geometry 

10,,) = 0.92 \d^P^P^PyPyd.) - \d^P.P.PyPyd,) ^ 3^^^ ^^^^3 

v2 

v2 

) = _0 Jo \diPxPxPyPy) + \PxPxPyPydj) 

-2\diPxPyPydj) + \diPxPyPydj) + \diPxPyPydj) 



+0.46 
-0.15 



\diPxPxPydj) - \diPxPxPydj) 



V2 

\d,^d,^n„n..fi..) 

+0.07 



\didiPxPyPy) - \PxPyPydjdj) 



V2 

_Q Q^ MiC^jPa^p'o^Py) - \PxPxPydjdj) 

+small terms . . . 



, \ , n on \diPxPxPyPy) - \PxPxPyPydj) 

>^o.) = +0.82 

\diPxPyPydj) - \diPxPyPydj) 



-0.31 



^/2 



_Q ryQ -'^\diPxPxPydj) + \diP^xPxPydj) + \diPxPxPydj) 

_Q ^ Jdid^PxPyPydj) + \diPxPyPydjdj) 

V2 



15 
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TABLE I. Exchange and hopping parameters in meV for a t — J Hamiltonian computed 

at different levels of c:alcnlation. 



Calculation level 


CASCI 


IDDCIl 


DDCI2 


Physics included 


Valence shell 


+dynamical 


+dynamical 






polarization 


polarization and correlation 


J' 


+5.48 


+4.14 


+4.87 


t' 


+37.5 


+28.7 


+44.2 




-1.05 


-4.64 


-5.04 


h 


-115.8 


-176.5 


-124.6 


J± 


-60.6 


-321.1 


-293.5 


t± 


-420.6 


-542.7 


-538.2 



TABLE IL Parameters in eV of the extended Hubbard Hamiltonian in the rung geom- 
etry. 



'■X 


ti 


,pd 
''X 


Ud - Vpd 


Up-S 


1.1 


2.0 


1.3 


3.6 


< .03 
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FIG. 1. Schematic structure of a NaViOf,, a) along the a-c plane, b) along the a-b 
plane. The oxygen atoms are are denoted by open circles, the vanadium atoms by filled 
circles and the sodium atoms by dots. 
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FIG. 3. Fragments of the crystal used for the cluster calculations, a) ty, and Jy ; b) 
t±, and J_i ; c) t' , and J'. 



21 




0.50 
0.25 
0.00 
0.25 
0.00 
0.25 
0.00 
0.25 
0.00 



" 1 

y 


1 

2a 


1 


_ ii 1 


2b 




_ . 1 


2c 


1 




26 

^ 1 


^ J 



-0.6 



-0.2 



0.2 



FIG. 4. Projected DOSS on the d magnetic atomic orbitals (1) and the p orbital of the 
bridging oxygon atom in the rung geometry (2). The crystal HF calculation (a). The rung 
cluster of this work (b). The rung cluster embedded with only one type of oxygen atom, 
all atoms supporting their formal charges (c). The rung cluster without embedding (d). 
The dotted vertical line shows the Fermi level. 
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FIG. 5. Equivalent magnetic system. The ellipsoids represent the delocalizcd magnetic 
electron on the rungs. The straight lines represent the ferromagnetic remaining interactions 
{J'eff) while the zig-zag lines represent the anti- ferromagnetic ones ( Jf^)- 
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